function dydt=AVP_ode_func(t,y,G,P,k)
u=-P.K*y;
dydt=P.A*y+P.B*u;
